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Abstract 

We studied screening masses in tlie equilibrium tliermodynamics of 
SU (2) and SU (3) pure gauge tlieories on tlie lattice. At a temperature 
of 2Tc we found strong evidence for dimensional reduction in the non- 
perturbative spectrum of screening masses. Mass ratios in the high 
temperature SU{2>) theory are consistent with those in the pure gauge 
theory in three dimensions. At the first order SU{3) phase transition 
we report the first measurement of the true scalar screening mass. 
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1 Introduction 



The equilibrium thermodynamics of a gauge theory is studied in the non- 
perturbative domain by lattice simulations of the partition function. Much 
is now known about the phase transitions in SU{2) and SU{3) pure gauge 
theory and in QCD, including the order, the transition temperature, Tc, 
entropy density, pressure, specific and latent heats and other such quantities 

i- 

Also of interest are the screening masses at finite temperatures. These 
are defined in general by the exponential spatial falloff of the correlation of 
two static sources. A classification of all masses is provided by the transfor- 
mation properties of the sources. For glueball-like screening masses, such a 
classification was performed in where the first measurements of several of 
these masses were reported. Extensive lattice data is available on screening 
masses from meson and baryon-like sources and on the correlation of 
Polyakov lines. 

A gauge invariant transfer matrix formulation is easy to write for the 
lattice regularised theory. For thermal physics, it is convenient to think of 
the transfer matrix in a spatial direction. The free energy and other bulk 
thermodynamic quantities involve only the largest eigenvalue of the trans- 
fer matrix. Correlation functions and the screening masses involve the ratio 
of this largest eigenvalue with specific other eigenvalues depending on the 
transformation properties of the source. Thus, the full spectrum of screening 
masses contains much more information about the theory than bulk thermo- 
dynamics can provide. 

A quantity of special interest is the Debye screening mass, Mjj, whose 
inverse gives the screening length of static choromo-electric fields. In the 
quantum theory, this screening mass plays an important role in regulating 
some infra-red singularities. Lattice measurements have, in the past, con- 
centrated on measuring a mass from the correlation of Polyakov loops, Mp. 
It was expected that when the gauge coupling becomes small, Mp = 2Md. 
Many years ago Nadkarni showed that the relation between Md and Mp is 
far from being so simple Subsequent lattice work focussed on methods 
of computing M^. 

Reisz and collaborators ||^, Q wrote down a dimensionally reduced theory 
which could be used to define Mo non-perturbatively. A recent paper [0 
used the representations of the symmetries of the transfer matrix to write 
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down operators whose correlations could be used to measure the Debye massQ 
and gave a general parametrisation of the perturbative and non-perturbative 
terms for M^. These parameters have since been determined in a lattice 
measurement of the Debye screening mass using a dimensionally reduced 
theory at very high temperatures [§]. 

A similar screening mass is required for the magnetic sector of the non- 
abelian gauge theory in the plasma, in order to get sensible results in the 
infrared. However, the magnetic mass is entirely non-perturbative in nature. 
It has been the object of many lattice studies [l^. In fact, a recent work 
[ p!0[| tries to extract Mr, and Mm from gauge fixed gluon propagators at finite 



temperature. 

Many years ago Linde discovered [O] that the thermal perturbation ex- 
pansion in non-Abelian gauge theories breaks down because the magnetic 
sector is not amenable to perturbative studies. A recent attempt to under- 
stand Linde's problem in the region where the coupling (? ^ 1 has invoked a 
sequence of dimensionally reduced effective theories [|12|. At length scales of 
1/ gT dimensional reduction yields a three dimensional SU{N) gauge theory 
coupled to an adjoint scalar field of mass Md ^ 0{gT). At longer scales, 
1/g'^T, the scalar field can be integrated out, and the leading terms in the 
effective theory correspond to a pure gauge theory in three dimensions. On 
the basis of this reduction it has been argued [0, |^ that a non-vanishing 
pole in magnetic gluon propagators is absent, and Linde's problem is finessed 
by confinement in the three dimensional gauge theory. 

In fact, dimensional reduction has often been used to explore finite tem- 
perature theories and it has long been argued that the dimensionally 



reduced theory is confining. The spatial string tension is known to be non- 
vanishing for T > Tc and scales as |T^. This, and other, lattice mea- 
surements reveal that the gauge coupling close to Tc are too large to trust 
perturbation theory. 

At temperatures of a few T^, the coupling g > I, and the length scales 
1/T ^ 1/gT ^ 1/gT'^, and hence cannot be decoupled. The construction 
of is not applicable. However, the transfer matrix argument assures us 
that Md is only the smallest in a hierarchy of screening masses. We can then 
use the spectrum of the screening masses to explore the symmetries of the 
transfer matrix and check whether or not dimensional reduction works; and 



^The group theoretical identification of Mp and A/^i was, in fact, first given in 
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if it does, then what is the nature of the reduced theory. 

Our major result is that the degeneracies of the spectrum of screening 
masses imphes that the symmetry group of the spatial transfer matrix is that 
of two dimensional rotations — implying dimensional reduction for T ^ 2Tc. 
The measured non-perturbative spectrum is similar to that in 3-d pure gauge 



theory 



The symmetries of the transfer matrix, and the physical consequences 
are discussed in Section 0. The group theory presented in this section is 
central to the rest of this paper. The extraction of screening masses in T > 
four- dimensional SU (3) and SU (2) pure gauge theories take up the next two 
sections. These may be skipped by those readers who are not interested in 
the details of the lattice simulations. Section |^ contains a summary of our 
lattice results and presents our conclusions on the nature of the dimensionally 
reduced theory. Several technical details are relegated to appendices. Ap- 
pendix gives the loop operators used in our computations. In Appendix 
^ we examine the possibility of explaining the mass spectrum in terms of 
perturbative multigluon states. Noise reduction techniques for the lattice 
simulations and the algorithm for projecting on to the lowest state in every 
channel are described in Appendices |C| and |^ respectively. 



2 Symmetries of the Transfer Matrix 
2.1 Group chains 

For the T = continuum Euclidean theory the symmetry of the transfer 
matrix is the direct product of the full rotation group 0(3) and the Z2 groups 
generated by charge conjugation, C, and time reversal. Irreps of 0(3) are 
labelled by the angular momentum and parity, J^, and of the full symmetry 
group by J^^ . For the lattice regularised theory, 0(3) breaks to the discrete 
subgroup of the symmetries of the cube. Oh. The consequent reduction of 



the irreps of 0(3) is well-known [T^ 



In the Euclidean formulation of the (continuum) equilibrium T > the- 
ory, the transfer matrix in one of the spatial directions is invariant under 
symmetries of the orthogonal slice. Such slices are three dimensional — two of 
which are spatial and one is the Euclidean time. The symmetry group is that 
of a cylinder, C = 0(2) x Z2. This Z2 factor is generated by cr^ : t i— » —t. 
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The non-Abelian group 0(2) contains the Abehan subgroup of rotations, 
S'0(2), and a 2-d parity, 11 : {x,y) i— > (x, —y). The 3-d parity P = C{7i)az, 
where C{n) is the rotation by vr in SO (2). 

In the high temperature hmit, the effective theory is expected to undergo 
dimensional reduction to a 3-d gauge theory. In such a theory, the transfer 
matrix must have the 0(2) symmetry of a 2-d slice. 0(2) has two one- 
dimensional irreps 0± (the 2-d scalar and pseudo-scalar) and an infinite tower 
of two-dimensional irreps M. Under the reduction of 5*0(3) to 0(2), the spin 
J irrep breaks as — 



J ^ On(j) + E ^here U{J) = (-1) 



(2.1) 



M=l 



On the lattice the irreps of C break further into irreps of the automor- 
phism group of a 2;-shce, the tetragonal group D\ = x Z2{P). For the 
three dimensional lattice theory, the automorphism group of the 2-d slice is 
C^, which is isomorphic to D^. A recent work on 3-d glueballs used this 
classification of the states [|15 |. 

This pattern of symmetry breaking is summarised by 



0(3) = S0{3) X Z2{P) 



C = 0(2) X Z^ia,) 



Oh = Ox Z^iP) 



X Z2{P) 



(2.2) 



0(2) 



Four of the Z2 factor groups are identical, and generated by the 3-d parity 
P. 0(2) and contain the 2-d parity 11. 



2.2 Point group representations 

In this paper we use the notation of for the crystallographic point groups 
and their irreducible representations (irreps). This subsection contains a 
discussion of the irreducible representations of the lattice symmetries. 

The group is generated by eight elements in five conjugacy classes — 
the identity {E), rotations of ±7r/2 around the 2;-axis (O4), rotations of vr 
around the z-axis (C*|), around the x or y-axes (O2) and around the directions 
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X ± y (6*2 = C4C2). The five irreps are called Ai, A2, Bi, B2 and E. The 
first four are one- dimensional and the last two-dimensional. The 10 irreps of 
Dj^ are obtained by adjoining the character of P to the irreps of D^. Under 
the reduction of Oh to D^, the irreps break up as — 



— s> Af, A2 — * -Bf, 

^ Af + 5f. (2.3) 

The notation Ai^2 is potentially confusing since these irreps can belong to 
both Oh and Df^. A similar confusion can also be caused by the notation 
which stands for the two-dimensional irrep of any group. In the rest of this 
paper we will mention the group involved whenever we mention one of these 
irreps. 

In the dimensionally reduced theory, the lattice symmetry group is ~ 
D4. Instead of calling the irreps by the names of the irreps, we denote 
them by the symbols A"*", A~ , B^, B~ and E respectively. The first four are 
one-dimensional irreps, and the sign which indexes them is the character of 
the 2-d parity 11 in the irrep. The breaking of irreps to is as — 



BtiB^ B^, B^,B2 B~, 



E^ E. (2.4) 

Our nomenclature for the irreps of is designed to be a mnemonic for 
the breaking of the Oh irreps (eq. p.3|) , but not for the breaking of to 

(eq. p.4| ). On the other hand, the naming convention followed in 0, |^ 
simplifies the latter at the expense of the former, since it proceeds from the 
isomorphism ~ x ^2(172). To make contact with the notation used 
elsewhere, note that P in [0] corresponds to our 2-d parity 11 and Rz to the 
character of the operator cr^, which is called P in 0. 

The 0(2) irreps break to as follows — 

0+ A+, and 0_ ^ A', 

(E (M = odd), 

M I B+ + B- (M = 2 mod 4), (2.5) 

[ + A- (M = mod 4). 

The A'^ corresponds to the scalar and the A' to the pseudo-scalar irrep of 
0(2). 
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Figure 1: Possible temperature dependence of the lowest screening masses in 
some channels. Dimensional reduction at high temperatures is signalled by 
a pattern of approximate degeneracies characteristic of symmetry. 

States or operators also carry a label C = ±1 for the charge conjugation 
symmetry. The charge conjugation operator reverses the direction of traver- 
sal of a loop, and hence takes the trace into its complex conjugate. C = 1 
states correspond to real parts of loops and C = —1 states to imaginary 
parts. For gauge groups with only real representations, such as SU{2), there 
are no C = — 1 states. The construction of irreps of from loops is given 
in Appendix ^. 

In finite temperature lattice simulations, the symmetry group of the trans- 
fer matrix is always D^. However, at low temperatures, we should expect to 
see an effective symmetry group Oh- At the other end of the temperature 
scale, if dimensional reduction is to be a good approximation, we should see 
an approximate symmetry. A guess at the temperature dependence of 
the screening masses is shown in Figure |l|. If the lattice is big enough, and 
the lattice spacing is sufficiently small, then we should see the more extended 
degeneracy of 0(2). In this case the A'^ and A~ of will be non-degenerate, 
since they correspond to different irreps of 0(2), but the and B~ will 
become degenerate, since they correspond to one irrep of 0(2). Differences 
between first and second order phase transitions should be visible in the Af~^ 
channel near T^.. 
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3 SU{3) Pure Gauge Theory 



We have simulated the SU (3) pure gauge theory with Wilson action at three 
temperatures. With A',- = 4 the critical coupling is Pc{Nr = 4) = 5.692 



[IS, |T9[. We have performed a simulation at T ^ with /? = 5.7. At 
T = 3Tc/2 the coupling is /?c(^t = 6) ~ 5.9 [|l^, and for T = 2Tc we use the 
coupling Pc{N^ = 8) fa 6.0 g. 

The simulations were performed with a Cabbibo-Marinari pseudo heat- 
bath update, where each update acted on three separate SU{2) subgroups 



by five Kennedy- Pendleton moves |21 . The class of loop operators mea- 



sured is listed in Appendix ^ Noise reduction involved a fuzzing procedure 
explained in Appendix y. For each irrep of the symmetry group, we pro- 
jected the measured correlation function to the ground state by a variational 
technique explained in Appendix 0. Successive measurements of correla- 
tion functions were separated by about one integrated auto-correlation time 
measured through the Polyakov loop. 

One of the techniques for extracting screening masses from correlation 
functions in the long direction (with sites) is to solve the equation 

cosh [m{z + l/2){Nj2-z-l)] ^ C{z + l) 
coAi[m{z + l/2){Nj2- z)] C{z) ^ ' ' 

for the local mass, m{z+l/2), given measurements of the correlation function 
C{z + 1) and C{z). The assumption that the correlation function is described 
by a single mass is borne out if there is a range of z for which the local mass 
is constant within errors. If there is such a plateau then we quote it as our 
estimate of the screening mass. 

In addition we have performed fits to correlation functions in the form 

^ i';;'"^^ + d - - --)' | . 0.2) 

[ cosh(yUo^z/2) cosh(/iiA2/2) J 

where r, /xq and /ii are fit parameters with the constraint /iq < /ii. Local 
mass estimates were deemed acceptable if the fitted value of /xq agreed with 
it. When the local masses were too noisy to show a plateau, we took /xq as our 
estimate of the screening mass. The goodness of the variational projection 
to the ground state was checked by observing whether the fitted parameter 
r was close to unity. 
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In minimising to perform the fits we took into account covariances of 
Since we normalised the variational correlator at separation zero to 



data 22 



unity, the error in this point is zero. However, the intrinsic variability in the 
zero distance correlator is then redistributed over all the points through the 
covariance. The maximum distance retained in the fit was always determined 
by the criterion that the correlation function at that distance should be more 
than 1-0" from zero. The estimates and errors of every measurable were 
constructed through a jack-knife procedure. 

3.1 T = Tc 



Operator 


(0,1) 


(0,2) 








r 


yUo 


x' 


r 


/"O 


AV 





5.7/6 


n 71 +0.04 


n 27+°-°i 

U-^ 1 -0.02 


6.3/6 


' "±-0.05 


27+0-01 

U-^ ' -0.02 


E++ 





2.6/6 


0.68t°:°^ 


n oQ+0.02 

U-'^<J-0.02 


2.0/6 


71+0-05 


n OO+0.02 
^•-^0-0.02 




T 


10.5/6 


0.811^^ 


0.26i°:°l 


10.9/6 


n oQ+0.03 


0.261°:°! 


Br 


T 






1.23 ±0.08 






1.2 ±0.1 



Table 1: Fitted masses in S'f/(3) theory for (3 = 5.7 on the 4 x 8^ x 16 lattice 
using 5000 configurations. Note that r depends on the basis set of operators, 
but /io does not. For the operator, a plateau in the local masses is used 
for the estimate of the screening mass. 

At /? = 5.7 we analysed configurations separated by 50 pseudo-heat-bath 
sweeps, discarding the first 20 configurations for thermalisation. Tunnelling 
between phases occurred every 600 sweeps on an average. The operators 
measured are listed in Appendix 0. Analyses were performed in subspaces 
corresponding to the Oh and D'^ irreps. 

From the results in Table |l|, it is clear that the lowest screening mass 
observed in the Ai~^ ^ and E~^^\q sectors coincides with that obtained in the 
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(b) 



-s — r 



(d) 



Figure 2: Af~^ correlation functions in SU{3) at P — 5.7 with (0, 1) variation, 
(a) Data and fits for the projection to the lowest mass (thermal Af~^: dark 
circles, T = A^^: grey circles, T = E'^'^: white circles), (b) local 
masses compared to the 1-a band around the best fit, (c) data and fits to 
the projection on the second mass, and (d) local masses compared to the 1-a 
error band on the fitted physical Af'^ mass. 



9 



Af^ sector. This is also clear from Figure 0. In addition, the eigenvec- 
tor of the variation over all operators was orthogonal (within errors) to the 
B^~^ ^ space and yielded the same mass as the A^^ ^. In the critical region, 

therefore, the spectrum of screening masses is organised in irreps of Df^. In 
[|] it was found that the screening masses could be organised into irreps of 
D\ at "iTcjl^ but not at Tc/2. Our observation extends this to the picture 
presented in Figure |l|. 

It is interesting to note that the projection to the A^^^ ground state, 
as measured by r, increases with the number of operators used. This is 
a generic feature of the variational method and clearly seen in Figure 0b. 
Another generic feature is visible in the same figure — the fitted mass is the 
same as the stable long-distance local mass. Also, because the fit uses all the 
data points, its error is slightly smaller than that of the local mass 

The lowest screening mass in the A\^ 
mates of the mass from Polyakov line corre 
mass", which goes to zero and decouples from the thermodynamics in the 
infinite volume limit. The mass in this channel, in the thermodynamic limit, 
should be finite, and on any finite lattice it would be the next-to-lowest A\'^ 
mass. We attempted to estimate this "physical mass" by solving the varia- 
tional problem for the second lowest eigenvalue. A common estimate of the 
mass is obtained from variation over all operators as well as only the A\ 

operators. This indicates that the "physical mass" in the A\ 
genuinely the lowest screening mass. Local masses and fits agree (see Figure 
0) and give a physical mass 

am(A++) = 0.86+[|;}J, (A^r = 4, /5 = 5.7). (3.3) 

Such a physical mass has also been estimated before for two 3-state spin 
models in three dimensions. An S'?7(3) spin model, obtainable from S'f/(3) 
gauge theory in the strong coupling regime, had physical mass 0.71 ± 0.06 



^ sector agrees with previous esti- 
ations ||2^. This is the "tunnelling 



T 

sector is 



2^ whereas the three state Potts model gave a mass of about 0.1 ||25|. The 



A{ mass measured at T = at /5c = 5.7 is 0.964 ± 0.012 ||2^, giving 

,2^ ' \ = 0.89 ±0.10. 3.4 

m(T = 0,A^-+) ^ ' 

In the scaling limit this ratio should be independent of /3c. 
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The Bi ^ operator is always a difference of loops. As a result the 
correlation function in this sector is much more noisy than in the Ai~^ sector. 
Nevertheless, we were able to follow the correlation function to distance 4, 
and obtain a plateau in the local masses. Our estimate of the screening mass, 
reported in Table 0, comes from the local mass, m(5/2) for both (0, 1) and 
(0, 2) variation. 



3.2 T = 3Tc/2 and 2Tc 

We made runs on 4 x 8^ x 16 lattices ai (3 = 5.9 and 6.1, corresponding 
to 3Tc/2 and 2Tc respectively. Since the integrated autocorrelation time for 
the Polyakov loop was estimated to be less than 10 sweeps, we analysed 
data separated by 10 sweeps after discarding the first 400 for thermalisation. 
5000 configurations were generated at 3Tc/2 and 10000 at 2Tc. At 2Tc we 
made two further runs. One was on a larger, 4 x 12^ x 16 lattice, where 
we collected 5000 configurations separated by 10 sweeps, after discarding 
the first 400 sweeps. The second was on a shorter 4 x 8^ x 12 lattice, with 
exactly the same statistics. We did not see any tunnelling events at all in any 
of the four runs; from a hot start the system quickly relaxed into one of the 
Z3 symmetric free-energy minima, and stayed there for the duration of the 
runs. Since measurements of autocorrelations of correlation functions showed 
that the integrated autocorrelation time did not exceed 1.5 measurements, 
these contributions to error estimates have been neglected in this section. 

We made measurements of the operators listed in Appendix 0. Each 
operator was replicated at five levels of fuzzing. In most channels we could 
follow the correlation functions out to distance 5 and found the (0, z) vari- 
ational ground state to be statistically well behaved even with z as large as 
3. The exceptions were the -Bf'*' and channels, which could be followed 
only to distance 3. As a result the variational ground state was stable only 
for z = 1 and 2. 

In each channel, the components of |0; z)t, the {0,z) variational ground 
state at temperature T, give the overlaps of the ground state with each 
operator. We normalised lO;^;)^ to unity in each jack-knife bin. We found 
that several components of this vector are numerically very stable from one 
jack-knife bin to another. The rest of the components fluctuate from bin to 
bin, and seem to fine tune the variational eigenvalue. The eigenvalue itself 
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is far more stable than any of the eigenvectors. 

We found that correlation functions at distance 1 differed qualitatively 
from the long distance correlation function in several respects, thus giving 
rise to certain systematics in the measurement of screening masses. 

• The overlap (0; 1|0; 3)2Tc differed significantly from unity. For example, 
this overlap was 0.73 ± 0.02 in the Af^ channel and 0.57 ± 0.05 for 
the ^2 On the other hand, the overlap (0; 2|0; 3)2t, = 0.90 ± 0.06, 
consistent with unity, for ^4^'*'. Similar results were obtained at 3Tc/2. 

• There was a strong effect on local masses. With (0, 1) variation, we 
usually found no plateau in the local masses. With (0, 2) or (0, 3) 
variations a plateau was often visible. 

• Fits to correlation functions also reflected this behaviour. No accept- 
able fit with one or two masses was found to the (0, 1) correlation 
function, whereas the (0, 3) correlator could be fitted with r ^ 1. 

In view of this, we quote results from the (0, 3) variational correlators in all 
channels except the -Bf^ and B2^, for which we quote results from (0,2) 
variation. All the screening masses we quote are obtained from local masses 
and verified by a fit. The exception is the ^2 channel which turns out to 
be noisy. In this case the quoted result is the fitted mass. 

The ground states in the Ai~^ , A^"^ and -B^"*" channels are temperature 
independent within errors. This is not so in the Af~^ and Bf~^ channels; the 
overlap of these two ground states at 3Tc/2 and 2Tc are both 0.58 ± 0.02. 
Interestingly, the masses are independent of the temperature. 

Our results for the screening masses are collected in Table ^. The most 
interesting result is the near equality of the Af~^ and ^2"^ screening masses 
at 2Tc. Similarly, the -Bf"*" and the -B^^ masses are equal and also equal to 
the Bj~^ mass. Perturbation theory cannot be used to explain this pattern of 
degeneracies because the P = 1 correlators require two gluon exchange, and 
the P = —1 correlation functions must have a minimum of four exchanged 
gluons. 

Finite volume effects are under good control, as shown by the three sep- 
arate runs on lattices of three sizes at 2Tc. The study in Q had shown that 
finite lattice spacing effects are also under control, by making two simulations 
at the same temperature but at two different lattice spacings. The values 
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0(2) 




3T,/2 


2T, 


irreps 


irreps 


small 


short 


small 


large 


o| 




0.64 ±0.01 


0.64 ±0.01 


0.65 ±0.01 


0.65 ±0.02 


o| 




0.80 ±0.02* 


0.61 ±0.07* 


0.69 ±0.04* 


0.68 ±0.04* 


01 




1.58 ±0.03 


1.61 ±0.06 


1.58 ±0.04 


1.56 ±0.04 


2+ 


5^ 


1.22 ±0.08 


1.22 ±0.05 


1.2 ±0.1 


1.27±0.06 


2+ 


5r+ 




1.2 ±0.1 


1.34 ±0.08 




2+ 




1.27 ±0.07 


1.3 ±0.1 


1.4 ±0.1 


1.3 ±0.1 



Table 2: SU{3) screening masses for T > Tc on small (4 x 8^ x 16), large 
(4x12^x16) and short (4 X 8^ X 12) lattices. The starred masses are obtained 
from two-mass fits, the rest from local masses. Data on correlators was 
not collected in two cases. 



of m/T in the Af^ and Ai~^ channels were seen to be independent of the 
lattice spacing. We expect that this is true also of the screening masses in 
other channels, but would certainly welcome a direct measurement. 

A comparison with zero temperature results is simple because T = 



measurements have been performed at both these couplings ^ ^ 
Using the results in 



m(3re/2. At 



we find that at P = 5.9 

m(3re/2,5++) 



0.78 ±0.05, 



m{T = 0,At+) 
Moreover, at /5 = 6.0 we find 

^ 1 ^ = 0.90 ±0.05, 



m(T = 0,E^ 
m(2T„ 5+^ 



0.88 ±0.09. (3.5) 



1.12 ±0.06. (3.6) 



m(T = 0,Al+) ■ ■ ' m{T = 0,E++) 

These comparisons are made between channels at T = which overlap the 
corresponding channel at T > 0. Although at this larger temperature the 
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T = Af^ mass is nearly equal to the thermal Af^ screening mass, the state 
is completely different. The eigenvector of the variational problem has equal 
overlaps with the T = Af^ and states. It is clear that the physics 
observed here is completely different from the T = physics. 



4 SU{2) Pure Gauge Theory 

The SU{2) pure gauge theory with Wilson action was simulated at three 
temperatures. With A^^ = 4 the critical coupling is ||2^ PdNr = 4) = 
2.2998. We performed two simulations close to — one with (3 = 2.30, and 
another at j3 = 2.25. With the lattice size we used, this other coupling is still 
within the critical region, as indicated by the Polyakov loop susceptibility 
[^. We also performed simulations at T = 2Tc with Pc{Nr = 8) = 2.51 and 
at T = 4Tc with p^iNr = 16) = 2.74 |0| . 



These simulations were performed with an over-relaxation |^ and a 



Kennedy- Pendleton heat-bath algorithm ||2T|. The class of loop operators 
measured is listed in Appendix |^. Each measurement of correlation func- 
tions was separated by about one integrated auto-correlation time measured 
through the Polyakov loop. The procedure for the analysis was identical to 
that for SU{3). 



4.1 T ^T, 



c 



For (3 = 2.30, we took one measurement on a 4 x 12^ x 16 lattice every 50 
heatbath sweeps, and worked with 10000 measurements after discarding the 
first 5000 sweeps for thermalisation. The set of operators measured is given 
in Appendix |^. Our results for masses are presented in Table |^ and the 
correlation functions and masses are displayed in Figure ^. 

The Af ^ masses coming from the Oh irreps Af and turn out to 
be identical. Note also that the parameter r is rather large, indicating a 
successful projection onto the ground state. The Af mass is significantly 
smaller than the T = mass at the same coupling, 1.22 ± 0.03 ||32[| . 

The Bf correlator is significantly more noisy. However the correlation 
could be followed out to distance five. The mass estimate is fairly stable, 
and the overlap with the ground state is rather good. This screening mass 
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1} i i 3 i 1, 6 1 H i i 3 i 1, 6 1 i 




Figure 3: Correlation functions in SU{2) at = 2.3 with (0, 1) variation, (a) 
Data and fits for the thermal correlator, (b) local masses compared 
to the 1-cr band around the best fit, (c) data and fits to the correlator, 
and (d) local masses compared to the 1-a error band on the fit. 
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/3 



2.25 



Operator 



At 
At 
At 
Bt 



T{A) 
T{E) 
T 

T{EX_ 



range 
[0:8 
[0:8 
[0:8 
[0:5 



0.32 
0.87 
1.08 
0.08 



0.48l°:°i 



0.54; 



-0.05 
-0.07 





U-J<J-0.05 



J^-^"-0.08 



range 
[0:7 
[0:8 
[0:8 
[0:5 



(0,2) 

0.11 
0.26 
1.16 
0.16 



Ato 



f^n+°-°^ 

U-OU_o.o6 



0.56: 



-0.06 
-0.06 



n ei4+0.05 
U-O"±-0.05 



1.11 



+0.10 
-0.10 



2.30 



At 
At 
At 
Bt 



T(A) 
T(E) 
T 



[0:8 
[0:8 
[0:8 
[0:5 



0.50 
0.63 
0.40 
0.50 





U-<J-L-0.01 

-0.01 
-0.01 



0.30 



O.SOt'of, 



-'-•-'-U_o.04 



[0:8 
[0:8 
[0:8 
[0:5 



1.02 
0.55 
0.36 
0.60 



0.32 



-0.01 
-0.02 



0.30: 



0.30: 



-0.01 
-0.01 



-0.01 
-0.01 



1.09t°:°l 



Table 3: Screening masses in the SU (2) theory near T^. 



is significantly higher than the At ^ screening mass, and much smaller than 
m(T = 0, E+) = 1.94 ± 0.08 [0, from which it comes. 

At P = 2.25 we took 10000 measurements separated by 50 sweeps after 
discarding the first 5000 sweeps. Our fits indicate that the two screening 
masses for At ^, those from the Oh irreps At and are almost degenerate 
at this point, although the common mass is larger than that measured at 
f3 = 2.30. The Bt ^ correlator was noisier, but good enough to yield a 
dependable measurement of the mass. This screening mass is clearly split 
from the At ^ mass coming from the same Oh irrep . Our results are 
summarised in Table ^ and are in accord with the expectation in Figure |l|. 
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4.2 T = 2Tc and iT^ 



We performed runs with 4 x 8^ x 16 and 4 x 12^ x 16 lattices at /3 = 2.51 
(2Tc) and /3 = 2.74 (4Tc). Three to five over- relaxation sweeps were followed 
by one heat-bath sweep. Since autocorrelation times of Polyakov loops and 
plaquettes were seen to be less than three such composite steps, we performed 
a measurement on every fifth composite step. In all these runs, the system 
quickly relaxed into one of the Z2 symmetric minima of the free energy and 
stayed there through the duration of the run. 

As explained in Appendix 0, we measured two sets of operators. The 
larger set. A, included all the operators in the smaller set, B. The marginal 
improvement in the measurement of screening masses did not compensate 
for the longer CPU time spent in constructing the extra operators. On the 
smaller lattice we measured only the set A. At 2Tc we made 10000 measure- 
ments on the small lattice and 10000 measurements of the operator set B on 
the larger lattice. At 4Tc we took 5000 measurements on the small lattice, 
20000 of the operator set B and 10000 of the operator set A on the large 
lattice. 

In the best cases we could follow the correlation function out to distance 
5, and the local mass m(3/2) already belonged to a plateau. The and 
B2 channels were noisy, and we could follow the correlator only to distance 
3. For these two channels as well, we quote m(3/2) as our estimate of the 
local mass. The A2 channel was more noisy than for the SU{3) theory. We 
were unable to make any measurement in this channel. Unlike the behaviour 
noticed for SU (3), the analysis from (0, 1) variation gave results in agreement 
with (0, 2) variation. Hence we report our analysis from (0, 1) variation. 

The eigenvectors corresponding to the ground states are very stable in ev- 
ery channel other than the B2, where it showed large bin-to-bin fiuctuations. 
The eigenvalues were always more stable than the eigenvectors. The ground 
state obtained by the variational method changed little between 2Tc and 4Tc 
for the Af , Ai and states; the eigenvectors had an overlap consistent 
with unity. However, in the Bi channel, this overlap was 0.70 ± 0.04. The 
screening mass m{Bj~) changes the most between 2Tc and 4Tc. 

Our results for the screening masses are summarized in Table ^. The 
measurement of the screening mass in the Af sector gives 




2.8 ±0.2 (T = 2TJ, 

2.9 ±0.2 (T = 4TJ. 



(4.1) 
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2T, 


4T, 


irrep 


small A 


large B 


small A 


large A 


large B 


At 


0.71 ±0.05 


0.69 ±0.04 


0.74 ±0.07 


0.73 ±0.05 


0.73 ±0.05 




1.14 ±0.02 


1.02 ±0.02 


1.10 ±0.05 


0.87 ±0.04 


0.88 ±0.04 


Bt 


1.18 ±0.03 


1.62 ±0.06 


0.69 ±0.02* 


1.08 ±0.06 


1.10 ±0.06 


Bi 


1.9 ±0.2 


1.8 ±0.1 


1.8 ±0.1 


1.55 ±0.06 


1.54 ±0.05 


Bt 


1.8 ±0.1 


1.75 ±0.05 


1.5 ±0.2 


1.60 ±0.08 


1.56 ±0.04 



Table 4: Estimates of screening masses in the SU (2) pure gauge theory from 
small (4 X 8^ X 16) and large (4 x 12^ x 16) lattices. The estimates are obtained 
from local masses, except the ones with a star. The latter are obtained from 
a fit. 



This is consistent with previous measurements from correlations of Polyakov 



loops, reported in The near-equality of the and the screening 

masses are characteristic of a dimensionally reduced theory, and cannot be 
explained in perturbation theory. 

The major difference between the SU(2) and SU{3) theories seems to 
be in the volume dependence of various screening masses. The Af mass 
is the only one which seems to be volume independent. The Bf screening 
mass increases with volume, and the rest decrease. A more extensive study 
is required to obtain the infinite volume limit of these masses. Only after 
this is done can we say more about the nature of the dimensionally reduced 
theory. 

Our measurements of m{T, Af) can be compared to the T = glueball 
masses measured at /? = 2.5 and 2.6 [^. Interpolation gives m(T = 0, Af) = 
0.70 ± 0.03 for p = 2.51. The mass ratio 



m(2T„A+) 
m(T = 0,At 



0.99 ±0.07, 



(4.2) 
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is close to unity. Thermal effects are seen in the ground state — the Af 
eigenvector has a large projection on both the T = Af and channels. 
Glueball masses at /3 = 2.74 for T = can be obtained by interpolating 
between the measurements at /3 = 2.7 and P = 2.85 |3^. We estimate 
m(T = 0, Af) = 0.35 ± 0.03. Then 



m(4T„ At) 
m{T = 0,At 



2.1 ±0.2, 



(4.3) 



indicating a large thermal shift. This shows that the physics of screening 
masses is quite different from that of T = glueball masses. 



5 Summary 

In this section we summarise the results of our lattice measurements and 
discuss the physics implied by it. We try to deduce some general features 
of the dimensionally reduced theory by comparing these results with what is 
known of three dimensional gauge theories. 

First we gather together our conclusions. We have found that for T > 
the spectrum of screening masses is completely consistent with a cylindrical 
symmetry of the spatial transfer matrix {D'l on the lattice). At temperatures 
of Tc and above, there is no remnant of the T = 0(3) rotational symmetry. 
The scalar representation of the cylinder group {Af^ on the lattice) gives the 
lowest screening mass at all temperatures in both SU{2) and SU{3) theories. 

Near Tc this lowest screening mass is very small. Since the SU{2) gauge 
theory undergoes a second order deconfining transition, it is expected that 
this screening mass should be precisely zero on infinite volume systems. The 
small non-zero value we observe can be ascribed to finite-size effects. 

The SU{3) gauge theory has a first order deconfining transition. Finite- 
sized systems near a first order phase transition show a small screening mass, 
which vanishes in the infinite volume limit, and is related to tunnelings, and 



hence the surface tension PBI, between the phases which coexist at a first 



order transition. Previous measurements had seen only this mass. A 



technical point of interest to specialists is that our measurement of this mass 
yields results consistent with previous observations. 

The important quantity for the study of screening masses is not this, but 
the finite screening mass — the "physical mass" in the Af^ channel for SU (3) 
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gauge theory. We have estimated it near Tc for the first time. 



T 


ot 




2+ 


2+ 


T 


3.4(4) 




4.9(3) 




3rp 

2 c 


2.56(4) 


6.3(1) 


4.9(3) 


5.1(3) 


2T, 


2.60(4) 


6.3(2) 


4.8(1) 


5.6(4) 



Table 5: The ratio m/T in various channels as a function of T in the SU{3) 
gauge theory. The Ai~^ screening mass shown near Tc is the "physical mass" . 
The masses are labelled by the name of the 4-d lattice irrep (upper line) and 
the 3-d continuum irrep (lower line). 

Our results for the screening masses are summarised in Table H. As dis- 
cussed in appendix the spectrum is not amenable to an understanding in 
terms of perturbative multigluon states, and gives nonperturbative informa- 
tion about the underlying effective theory. 

Evidence for dimensional reduction at 2Tc comes from degeneracies in the 
spectrum of the transfer matrix. In the SU{3) theory m{Af~^) ~ m{Ai~^), 
and m{Bi^) ~ m{B2^) ~ m{Bi^). The first two sets of equalities are suffi- 
cient to argue for dimensional reduction of the lattice cutoff theory. The third 
equality implies 0(2) symmetry and hence allows us to make the stronger 
statement that the lattice artifacts are small. In the SU{2) theory, although 
we have not eliminated lattice artifacts, dimensional reduction is shown by 
the relation m(i?^) ~ m{Bi). 

Now we turn to possible interpretations of our detailed observations. The 
numerical values of the screening masses we have observed constrain the 
form of the three dimensional effective theory that describes equilibrium 4- 
d thermal gauge theories. In the scaling region of 3-d SU{N) pure gauge 
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theories, the glueball mass ratios 



. 0.60, ^ 0.78. (5.1) 

are almost independent of for gauge groups SU{N) Our measure- 

ments at 2Tc in the SU (3) theory give 

m(A++) m(5++) , , 

— ^— -( = 0.54 ±0.02, — T-r-TV = 0.76 ± 0.02. 5.2 

m{B++) m{A-+) ^ ' 

From the spectrum it seems hkely that the dimensionally reduced theory 
corresponding to the finite temperature SU (3) theory may be a 3-d pure 
gauge theory. 

One final observation — the A\ mass in the SU{2) theory agrees with 
many other measurements (performed through Polyakov loop correlations), 
and is expected to be independent of the lattice spacing. It is also seen to be 
independent of the lattice volume. In addition, it agrees numerically with the 
Af^ screening mass observed in the SU{3) theory. It would be interesting 
to study whether the infinite volume limit of the other screening masses in 
both these theories show a similar agreement. 

We will report on several technical points in future. A study of finite 
lattice spacing effects is under way. We are also performing a more extensive 
study of finite volume effects. We have not studied the two-dimensional 
irreps, E, of D^, but expect that the screening masses in these channels 
would add to our understanding of this problem. 

We would like to thank Rajiv Gavai for discussions. 
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A Loop Operators 



We specify a loop with the notation Ui{fi, u, X,- ■ ■). This denotes a product 
of hnk matrices starting with f/^(i) and proceeding along the links in the 
directions u, A, etc. The loops used in this work are drawn from the set — 

= Re Tr f/j(/i, z/, — /i, — i^), 

= Re Tr f/i(/i, /i, z/, — /i, — /i, — z/), 

= Re Tr f/j(/i, u, p, -/i, -u, -p), 

= ReTr Ui{p,u,-p,p,-u,-p), 

Ol = Re Tr ?7i(/i,/i, z/,p, -/i, -/i, -p, -z/), 

Of = Re Tr Ui{p, p, z/, -p, u, -p, -z/, -z/), 

O4 = ReTr Ui{p,p,u,-p,p,-p,-p,-u), (A.l) 

05= ReTr Ui{p,p,iy,-p,p,-i^,-~p,-p), 

O7 = ReTr Ui{p,p,u,p,-u,-p,-p,-p), 

= Re Tr f/i(/i, z/, -p, -z/, -/i,p,/i, -p), 

Ow = Re Tr Ui{p, u, -p, -u, -p, -z/, p, z/), 

= Re Tr Ui{p, z/, -p, -z/, -p, z/, p, -z/), 

= Re Tr f/,(p,p, z/, -p,p, -/i, -z/, -p). 

Our convention for naming these loops follows that of |T^. The plaquette 
(P^), 6-link planar (P^), twisted (T^) and bent (P^) loops, and the 8-link 
loop 0^4 were considered in detail earlier We have also used the double 
traversal of some of these loops — 

= ReTiUiUi, where O = Re Tr t/„ (A.2) 

where Ui is the SU{N) matrix corresponding to a loop. The representation 
content of such pairs O and are identical. 

The irreps of the symmetry group can be constructed by acting on any 
loop by the projection operators of the group. For we use the operators 

Af = {E + C^){E + Ci){E + C2){E±P), 

Af = {E + Ci){E + Cl){E-C2){E±P), 

Bf = {E-C^){E + Cl){E + C2){E±P), (A.3) 

Bf = {E ~ C^){E + Ci){E - C2){E ± P). 
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We have not used the remaining four projectors, which are the two indepen- 
dent sets each of and E~ projectors that can be obtained by changing 
the sign of C| in the above formulae. The irrep content of the loops in eq. 
( [A. 1|) can simply be obtained using these projectors in a small Mathematica 



program. 

In the SU (3) measurements we used the following set of operators- 





(2) 


P^ (2) (P4)2, (3) P\ 


(3) {pr, 


(1) 0?4, (1) Of,, (1) 0?„ 




(3) 


ofs, (3) {o%y 








(1) 


01, (1) Of,, (3) Of,, 


(3) {of,r 






(1) 


B\ (1) {B^f, (1) Of,, (1) Of,, ( 


1) 0?„ (3) Of,, (3) (0?3) 




(1) 


p^ (1) {p'f, (3) p\ 


(3) {py, 


(1) Of,, (1) 0?e, 




(3) 


Ofs, (3) {Of,r 








(1) 


of,, (1) Of,, (3) Of,, 


(3) {of,r 






(1) 


Of,, (1) Of,, (3) Of,, 


(3) {of,r 





For the simulation near Tc only the plaquette and 6-link planar loops were 
used. 
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For the SU (2) measurements a bigger set of operators was used- 



At 

Al 
At 
A2 

Bt 



Bt 
B, 



(2) p^ 

(2) 01 

(3) 01 
(3) 01 
(1) B\ 
(1) 0?6, 
(1) P^ 
(3) 01 
(3) Of, 
(1) 5^ 



3) (2) 5^ (1) T^ (2) (2) o|, (3) 01 

3) O?, (2) 01, (2) (2) (3) Of, 
2) 0«, (2) 01, (2) 0?6, (3) 0?3 



1) 01, (1) O?, (1) 01, (1) 0?e,(3) O 
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1) 01, (1) Of, (3) 01 (1) Of, (2) 0«, (1) 01 



14' 



(1) (3) O 
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(3) (1) P^ (1) Of, (1) Of, (3) 01 (1) Of, 
;i) 01, (1) 0?„ (1) Of,, (3) 0?3 

;i) o«, (1) 01, (1) o?„ (3) Of, 

;i) T^ (1) Oi, (1) Of, (3) 01 (2) 01, (1) o?, 
(2) 0?4, (2) Of,, (1) (3) Of, 

(1) P^ (1) (1) Of, (3) 01 (2) 01, (2) O?, (2) 0?4, 

(2) Of,, (1) (3) Of, 

along with the double traversals of each of them. This full set is the SET A 



of Section IJ. Set B contained only the operators P^ P^ P^ T^, Of, Of^, 



of, and Ofg. For the simulations near Tc only the plaquette and the six-link 
loops were used. 

B Multi-gluon States 

The simplest way of understanding the screening masses in the high temper- 
ature phase would be to accomodate it in a perturbative framework. Corre- 
lators of loops will then be dominated by suitable multigluon states, and the 
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screening masses would give information on Md and Mm 0- Fortunately, 
such a hypothesis is open to a direct numerical test. In this appendix we show 
that the screening masses cannot be understood in perturbation theory. 



k 


At 


A A 


(0,0,0) 
(0,k,0) 
(k,0,0) 
(0,k,k) 
(0,k,k') 
(k,k',0) 

(k,k',k") 


At, 

A^,Bi,E+ 
A^,A,,B^,B^,2E+ 
At,A,,B+,B,,E+,E- 
At,At,Bt,Bt,E^,E^ 


E~ 

At,Bt,E- 
E+,E- 

At,At,Bt,B+,2E- 
At, A2, Bi, B2, 2E^ 
At,A2,Bt,B2,E+,E- 
Af,At,Bt,Bt,E^,E^ 



Table 6: The representation content of gluon operators at various points, 
k, of the Brillouin zone. The first component of k is /cq- Note that all these 
irreps come with C = —1. 

We first classify the irreps of Df^ which can be obtained by specific multi- 
gluon exchanges. The gluons carry arbitrary momenta, and are classified by 
representations of the space group. We break these irreps of the space group 
under the point group. Then the allowed representations for colour singlet 
zero momentum correlators are obtained by combining the gluon represen- 
tations by the Clebsch-Gordan series for the point group and applying 
appropriate exchange symmetries. 

We begin by specifying the lattice analogue of gluon field operators in 
momentum space by the Fourier transform of a projection of link matrices 
onto the SU (N) algebra — 

^,(k) = ^^e'''-[f/,(x)-f/t(x)-ImTrf/,(x)]. (B.l) 

X 
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Here x takes values in a slice orthogonal to the z-direction, and k in the 
corresponding Brillouin zone. In general, the action of carries one k into 
another. The representations built over the orbit of k under the action of 
D'^ are usually large and reducible. The representation content of such gluon 
field operators is gauge independent, and is shown in Table 



k 


^(2) 


r^(2) ^(2) 


(7(2) G(2) 

XX , yy 


^(2) 
xy 


(0,0,0) 


At 




E+ 


At,Bt 


Bt 






(0,k,0) 


At, 


Bt 


E+ 


At,Bt 


4 + 

^2 ) 


Bt 




(k,0,0) 


At 




E+ 


At,Bt 


Bt 






(0,k,k) 


A + 
^1 ) 


Bt 


2E+ 


4+ 4 + 

^1 5 ^2 5 


At, 


Bt 












Bt,Bt 








(0,k,k') 


At, 


At, 


2E+ 


A+ 4 + 

^1 5 ^2 5 


At, 


At 






Bt, 


Bt 




Bt, Bt 


Bt, 


Bt 




(k,k',0) 


A + 

^1 7 


Bt,E+ 


At,Bt,E+ 


At,Bt,E+ 


4 + 

^2 ) 


Bt 


,E+ 


(k,k',k") 


At, 


A+ R+ 


A+ A+ f>+ 


At , At 1 Bt 1 


At, 


At 


Bt, 




Bt, 


E+ 


Bt,E+ 


Bt,E+ 


Bt, 


E+ 





Table 7: The irreps for symmetric colour singlet C = 1 two-gluon states with 
total momentum zero are listed for the momentum, k, of one of the gluon at 
various points in the Brillouin zone. 

Gauge invariant C = 1 states of zero momentum can be constructed as 
linear combinations of the composite operators 

G^i) = Tr[A,{k)A^i-k)], (B.2) 

where the trace is over SU{N) generators. Correlators of such an operator 
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will decay, in leading order, with the mass E^(k) + -E^(k), where Eq = Ed, 
Ei,2 = Em and 

2 

smh\EDM/2) = smh\MD,m/2) + J2sm\ki/2). (B.3) 

i=0 

The cyclic property of traces ensures that G*^^-* is symmetric under any opera- 
tion that flips the polarisation indices and simultaneously changes the sign of 
k. In addition, if we require that the state be symmetric under the exchange 
of the gluon fields, then only P = 1 irreps are allowed. The representation 
content of these operators in all parts of the Brillouin zone is given in Table 

0- 

In order to obtain P = — 1, C = 1 irreps we have to go to combinations 
of four gluon field operators. Colour singlet gauge invariant correlators with 
C = —1 start with composite operators of three gluon fields. 

The reduction of multi-gluon operators then tells us that 

• For C = 1, the lowest mass observed in any P = 1 channel would be 
at most half of the lowest mass observed in any P = — 1 channel, since 
the former are obtained by two-gluon exchange but the latter require 
at least four-gluon exchange. 

• If Em{0) < Ed{0), then m(A++) = m(P++) = m{B++). 

• If Ed{0) < Em{0) < Ed(K), for K = {0,27r/N^,0), then m(Pi++) = 
m{B2~^) and the mass difference between the Af~^ and Bf^ states is 
2[Em{0) — -E'd(O)]. Here is the spatial size of the z-shce. 

• If -^^(0) > -Ed(K), then m{Bi~^) < m{B2^) and the mass difference 
mlBi'^) — m{Ai~^) decreases with lattice size. 

The first condition is clearly violated in the spectra we obtain for both 
SU{2) and SU{3) theories. At 2Tc in both the theories we found m(Pf"'") ~ 
m{B2~^). In addition, in the SU{3) theory m{Af~^) ~ m{A2~^). Both these 
observations violate the first condition. This is sufficient evidence for a fail- 
ure of the interpretation of these masses as screening masses for perturbative 
multigluon states. 

If we ignore this, and restrict ourselves to the P = 1, C = 1 states 
only, then one might be tempted to match the SU{3) spectrum to the third 
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condition above, since the mass difference m{Bi'^) — m{Ai^) is non-zero 
and independent of the lattice size. However, from eq. (|B.3|) , the condition 
Mm < -^^(k) is satisfied for A^^ = 8 but not for A^^, = 12. Thus, none of 
the conditions above hold, and we conclude that the spectrum of screening 
masses cannot be given a perturbative interpretation of multigluon states, 
and indeed gives information about the nonperturbative spectrum of the 



effective theory. A similar statement has been made in ref. |33] in the 
context of the massgap obtained from Polyakov loop correlators at moderate 
temperatures, on the basis of the detailed decay patterns of the correlators 
seen in lattice studies. 



C Improved Operators 

Loop correlations are known to be very noisy. In order to increase the sig- 
nal/background ratio, we used a hybrid of Teper's doubled-link fuzzing pro- 
cedure ||3^ and the smearing procedure adopted by the APE collaboration 
. We define fuzzed links at level / + 1 recursively in terms of those at level 
/ by the equation 

maxtr(Mt^f/^;;^)) where M,,, = U^l + uS^uj^^^^^:^,,, (C.l) 

jjH+i) g^j-Q elements of SU{N), and the links for / = are those generated by 
the Monte Carlo procedure. For general SU{N), this maximisation is most 
easily accomplished using the "polar" decomposition of a general complex 
matrix to write 

M = iuU'-^+^^H, (C.2) 

where u is a complex number of unit modulus, H is hermitian and [/('+^) is 
special unitary. There is a discrete ambiguity in this decomposition, corre- 
sponding to the signs of the eigenvalues of H. When all the eigenvalues of H 
are chosen to be positive, U^''^^^ maximises the trace in eq. For SU{2) 

the algorithm is simpler since if is a multiple of the identity. The projection 
then involves only a division of M by the square root of its determinant. 

In a test run with SU{3) at (3 = 5.7 on a 4'^ x 12 lattice, the procedure in 
eq. ( |C.1| ) was found to perform better than doubled-link fuzzing. Since the 

^We would like to thank Gautam Mandal and Avinash Dhar for a discussion of this 
point. 
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latter technique is known to work well on larger lattices, we conclude that 
the problem is due to the fact that with small lattices, only a small number of 
doubled-link fuzzing steps is possible. Presumably on larger lattices, where 
more fuzzing levels can be reached, equally good results can be obtained with 
either fuzzing technique. For the SU{3) theory we worked with eq. ( |C.1|) and 
/ < 4. 

For SU{2), since we use a lattice which has 12 spatial sites, upto three 
levels of doubled-link fuzzing can be performed for the spatial links. For finite 
temperature problems it is perfectly all right if the number of fuzzing steps 
in the time direction is different from that in other directions. We perform 
only one doubled-link fuzzing in the time direction. We checked that using 
three levels of doubled-link fuzzing gave a better projection than seven steps 
of (C.l), and therefore used the former, for our runs near T^. For our runs 
at 2Tc and 4Tc, we experimented with a combination of doubled-link fuzzing 
and (C.l); using a combination of one doubled-link fuzzing followed by two 
steps of (C.l), the whole set being repeated once. Since this gave a slight 
improvement over three steps of doubled-link fuzzing, we used this technique 
at these higher temperatures. For our runs on the smaller lattice we used 
one doubled-link fuzzing followed by 4 steps of eq. (|C.1|) . 

D Variational Correlators 

It is not known a priori which linear combination of loop operators acting on 
the vacuum generates the state with the lowest mass in a channel with given 
quantum numbers. However, such a state will give a correlation function 
which has the slowest possible decay with increasing separation. Given a 
basis set of loop operators, we can try to construct a linear combination 
which satisfies this property of slowest decay. This is the idea of a widely used 
variational technique [^]. Since we have found no discussion in the literature 
of a numerically stable algorithm for its implementation, we document such 
a method here. 

We construct cross correlations between all the loop operators at our dis- 
posal to yield the (symmetric) matrix of correlations Cij{z). A combination 
of operators which has large projection to the ground state is obtained by 
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solving the variational problem over Y 

'Y^C{zi)Y 



max 



Y^C{zo)Y^ 



\{Zq,Zi), {Zq<Zi). (D.l) 



If C{zq) is positive definite, as guaranteed by the reflection positivity of 
the Wilson action, then this extremisation problem reduces to finding the 
maximum eigenvalue of the system — 

C{z,)Y = \{z,,z^)C{zo)Y. (D.2) 

With the corresponding eigenvector, we define the (zq? zi) variational corre- 
lator 

C,,,,,{z) = Y^C{z)Y. (D.3) 

We can utilise the freedom of normalising the eigenvector Y to set the vari- 
ational correlator to unity at separation ^ = 0. 

Reflection positivity of the action guarantees that C{zo) is positive defi- 
nite. It can be treated as a metric, and after appropriate scaling, the problem 
in eq. ( D.l ) can be phrased as the extremisation of a quadratic form over a 



sphere — leading to the usual matrix eigenvalue problem ||4l|. Algorithmi- 



cally, this naive idea can be implemented by transforming both sides of eq. 
(p.2|) to the basis where C{zq) is diagonal, absorbing the diagonal elements 
into Y by appropriate rescaling, and then solving the usual eigenvalue prob- 
lem for this transformed C{zi). However, if some of the eigenvalues of C{zq) 
are small, then the extremum problem is ill-conditioned because the solution 
is sent off to infinity along the nearly fiat directions. 

With finite statistics the problem may be even worse. Due to statistical 
fluctuations, the measured correlation matrix may not be positive definite. 
It is then better to treat the problematic directions as exactly fiat, since 
this discards the subset of the data which is most corrupted by noise. The 
solution is easily specified by going to the basis in which C{zq) is diagonal 
and blocking the matrices into the form 



where the eigenvalues, Cj, of C{zo) which satisfy the cut condition 

Ci < eeo (D.5) 
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have been set to zero. Here eo is the maximum eigenvalue of C{zo)- The 
diagonal sub-matrix G is positive definite and has a condition number less 
than 1/e. In this basis, eq. (p.l|) is equivalent to the set of equations — 



CnY, + CuY2 = XGY, and Cf^Y, + C22Y2 = 0. (D.6) 

Solving the latter for Y2 and substituting into the former gives the eigenvalue 
problem — 

Cii — C12C22C12 Yi = AGYi, (D.7) 

which is well-defined and numerically well-conditioned as long as C22 is in- 
vertible. 

Notice, however, that eq. (p.2|) is ill-conditioned only if both C{zi) and 
C{zo) have nearly fiat directions; otherwise their roles may be interchanged 
by inverting eq. ( p.l|) and converting the maximum problem into that of 



finding a minimum. If this cannot be done, then we must take care of the 
case that C22 is not invertible. 

The m X n matrix is a map from to (i.e., Yi has n real 
components and Y2 has m). Its range is the subset of to which the whole 
of -R" is mapped. If C22 is singular, then its null-space is contained in the 
complement of the range of C^. As a result, the second of eq. ( p.6|) can 
be solved by a left-multiplication by any matrix which coincides with the 
inverse of C22 in the complement of its null-space (i.e., in the range of C^). 
Therefore, in eq. (p. 7] ) C22 can be replaced by a pseudo-inverse 



C22^ = V^k-^V where C22 = V^AV, (D.8) 

V is an orthogonal matrix, A is diagonal, and the small components of A 
are set to zero in the pseudo-inverse A~^. With this definition of C22, the 
generalised eigenvalue problem is completely well-defined and eq. ( p.7| ) may 
be solved by the naive algorithm described earlier. 

A very nice numerical illustration is provided by our data on the 95 x 95 
matrix of the Af~^ correlation function at 2Tc for the SU{3) gauge theory 
(section P^). The distribution of 



^ = logio A , (D-9) 




is shown in Figure The ten small eigenvalues clustered at the end of a 
huge spectral gap have both positive and negative signs, and are due to noise 
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Figure 4: Histogram of normalised eigenvalues of the 95 x 95 correlation 
matrix at distance zero, C(0) in the Af^ channel at 2Tc for the SU{3) theory. 



in the data. Because they are so well separated, the cut e (in eq. p.5|) can be 
chosen to have any value between lO""^ and 10^^^. The results for eigenvalues 
and eigenvectors are stable in this whole range of choices. Somewhat smaller 
values, e ~ 10~^, are also found to be acceptable and are in fact preferred 
for reasons of numerical stability of the linear algebra routines. Removing 
the cut destabilises the problem completely. The spectral distributions are 
similar in most channels for both SU{2) and SU{3) theories. In a few cases 
the spectral distribution is gapless but has a long tail. Even in such cases, 
e ~ 10~^-10~^ give stable results. 
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